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Abstract. Damping of magnetic fields via ambipolar diffusion and decay of magnetohy- 
drodynamical (MHD) turbulence in the post decoupling era heats the intergalactic medium 
(IGM). Delayed recombination of hydrogen atoms in the IGM yields an optical depth to scat¬ 
tering of the cosmic microwave background (CMB). The optical depth generated at z 3> 10 
does not affect the “reionization bump” of the CMB polarization power spectrum at low 
multipoles, but affects the temperature and polarization power spectra at high multipoles. 
Writing the present-day energy density of fields smoothed over the damping scale at the 
decoupling epoch as ps,o = -Bq/2, we constrain Bo as a function of the spectral index, ub- 
Using the Planck 2013 likelihood code that uses the Planck temperature and lensing data to¬ 
gether with the WMAP 9-year polarization data, we find the 95% upper bounds of Bq < 0.63, 
0.39, and 0.18 nG for ns = —2.9, —2.5, and —1.5, respectively. For these spectral indices, 
the optical depth is dominated by dissipation of the decaying MHD turbulence that occurs 
shortly after the decoupling epoch. Our limits are stronger than the previous limits ignoring 
the effects of the fields on ionization history. Inverse Compton scattering of CMB photons 
off electrons in the heated IGM distorts the thermal spectrum of CMB. Our limits on Bo 
imply that the y-type distortion from dissipation of fields in the post decoupling era should 
be smaller than 10 “ 9 , 4 x 10 9 , and 10 9 , respectively. 
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1 Introduction 


Damping of magnetic fields affects the thermal and ionization history of the universe. The 

damping processes before the decoupling epoch can be approximately modeled by damping 
of three magnetic modes. Of these, damping of the fast magnetosonic waves proceeds in a 
similar way as the Silk damping of density perturbations, resulting in a damping scale on 
the order of the Silk scale. Slow magnetosonic and Alfven waves can avoid damping up to 
smaller scales, resulting in a larger damping wave number. The value at the decoupling epoch 
is given by [1, 2] 


kd,dec 


299.66 / B 0 \ _1 
cos 9 \ 1 nG ) 


Mpc -1 , 


( 1 . 1 ) 


for the best-fit ACDM model of the “Planck 2013+WP” data [3]. Here, 9 is the angle between 
the wave vector and the field direction which will be set to zero in this paper. 

We assume that the magnetic held is a non-helical and Gaussian random held with its 
two-point function in Fourier space given by 


{B*(k)Bj(q)) = (2vr ) 3 S(k - q)P B (k ) , 


( 1 . 2 ) 


where the power spectrum, P B (k), is assumed to be a power law, P B {k) = A B k nB , with the 
amplitude, A B , and the spectral index, n B . This can be expressed in terms of the ensemble 
average of the present-day magnetic energy density by 

(PB,o) = j |^P B ,o(fc)e“ 1 2 k) , (1.3) 

where k c is a Gaussian smoothing scale. The power-law magnetic power spectrum then gives 


4*2 2 (»b+3)/ 2 (k ye 

= irypsiy (fc) 


(1.4) 


For convenience, we dehne a smoothed magnetic held strength at the present epoch, 
B\ c 0 i as (p B ,o) = \P\ C 0 ) where A c = 2-k/ k c is a wavenumber corresponding to k c . In this 
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paper, we shall set k c to the maximal wave number determined by the damping scale at the 
decoupling epoch, k^deci given by equation (1.1). For simplicity we shall write Bq = B\ ddect o 
throughout the rest of the paper. 

Primordial magnetic fields present before the decoupling epoch can generate scalar, vec¬ 
tor, and tensor modes of CMB temperature and polarization anisotropies [4-8]. In general, 
there are two different types; namely, the compensated magnetic mode which is similar to 
an isocurvature mode, and the passive mode due to the presence of anisotropic stress of the 
magnetic held before neutrino decoupling. The latter is just like a curvature mode with its 
amplitude determined by the magnetic anisotropic stress and the time of generation of the 
primordial magnetic held. These can be important contributions to the CMB anisotropies, 
depending on the magnetic held parameters. On small angular scales, the scalar compensated 
mode as well as the vector mode for nG helds with positive spectral indices, n B > 0, can dom¬ 
inate over the primary CMB anisotropies. On large angular scales, the tensor compensated 
mode for a few nG helds with a nearly scale invariant spectrum, ng = —2.9, can dominate 
[9]. However, here we shall ignore these contributions and derive constraints focusing only on 
the effect of dissipation of helds in the post decoupling era. Since we only consider n B < 0, 
the former contributions can be ignored. The limit we obtain for ng = —2.9 is smaller than 
nG, making the latter contributions sub-dominant. In any case, including these contributions 
in the analysis would make the constraints on Bq stronger. 

After the decoupling epoch, magnetic helds are damped by ambipolar diffusion and 
decaying MHD turbulence [10]. This leads to a change in the ionization and thermal history, 
hence the visibility function of CMB anisotropies. In ref. [11], we have calculated the thermal 
and ionization history of the IGM including dissipation of helds. Ionization gives an optical 
depth to Thomson scattering of the CMB. By writing the total optical depth as Ttot(Bo,n B ) = 
Ttot(Bo = 0) + Ar(Bo,ns), the additional contribution from dissipation of magnetic helds 
Ar is given by 1 

Ar(B 0 ,n B ) = 0.0241 (-n B )-° 0370 

x 5.2815xl0“ 12 (—n s ) 23 - 8731 +5.4xl0- 3 (^) 3 ' 3706 —7.1xl0" 3 (—ns) 1 - 948 ^) 2 ' 071 ^! ^ 


for ng < 0. 

In general, there are two important contributions to the total optical depth out to the 
decoupling epoch: the contribution due to reionization, r re j 0 , and that generated at high 
redshifts close to decoupling, Thi g h- Z - Hence T to t = Treio"^ T high-z- Whereas r re j 0 determines 
the CMB polarization angular power spectrum at low multipoles (“reionization bump”), the 
total optical depth out to decoupling determines the overall suppression of the temperature 
power spectrum at high multipoles as Cj T —>• Cf T e~ 2rtot . As the contribution to the optical 
depth from dissipation of helds comes from high redshifts, 2 3> 10, it has little effect on r re j 0 
and hence the CMB polarization power spectrum at low multipoles. Therefore, the change 
in the amplitude of the temperature angular power spectrum induced by the magnetic held 
dissipation is determined by At(Bq, ng), and the temperature data at high multipoles will be 

1 This formula is a minor update of the formula we gave in ref. [11]. We now use the Planck 13+WP 

best fit cosmological parameters. Following ref. [9, 12], we calculate the photoionization rates at the photon 
temperature instead of the electron temperature in the modified version of Recfast++, and include the colli- 

sional ionization. Finally, we have also corrected a numerical error in the amplitude of the ambipolar diffusion 
heating rate. We are grateful to J. Chluba for helping us implement these updates. 
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the key to measure this effect. In this paper, we shall use the Planck 2013 temperature data 
[13, 14] (including the lensing data [15]) and the Planck 13 polarization likelihood [14] derived 
from the WMAP 9-year polarization data [16, 17] to constrain B$ for some representative 
values of ns- In particular, the lensing data play an important role in isolating the optical 
depth in the temperature data. 

The rest of the paper is organized as follows. In section 2, we calculate the CMB 
temperature and polarization power spectra with a modified thermal and reionization history 
due to the dissipation of the magnetic field in the post decoupling era. In section 3, we obtain 
constraints on the magnetic field strength for some representative values of ns- We conclude 
in section 4. 


2 CMB anisotropies with dissipation of magnetic fields in the post-decoupling 
era 


Dissipation of magnetic fields in the post decoupling era takes place by two processes; namely, 
ambipolar diffusion and decaying MHD turbulence [10] . These processes affect the evolution 
of the electron temperature, T e , as 


8/OyCTT 


T e = —2—T e + " c (T y - T e ) + 

a 1 + x e 3m e c ' 7 


x P T 


1.5 ksfle ’ 


( 2 . 1 ) 


where the dissipation rate, T = Ti n + Tdecayi includes the contribution from ambipolar diffu¬ 
sion, Tin, and decaying MHD turbulence, Ld ecay . The other variables are: a the scale factor, 
T 7 the photon temperature, ax the Thomson scattering cross section, x e the ionization frac¬ 
tion, n e the electron number density, p 7 the photon energy density, m e the electron rest mass, 
c the speed of light, and ks the Boltzmann constant. 

Ambipolar diffusion is due to different velocities of the neutral and the ionized matter 
components in the presence of a magnetic field. As this velocity difference is caused by 
the action of the Lorentz force on the ionized component, the dissipation rate of ambipolar 
diffusion is determined by the Lorentz term as [10] 

T in = i6 J^ 2 ^ \(VxB)xB\ 2 , (2.2) 

where p n , pi , and pb are the energy densities of neutral hydrogen, ionized hydrogen, and the 
total baryons, respectively, and 7 is the coupling between the ionized and neutral component 
given by 7 ~ (ov) H + H /(2m#) with {crv) H + H = 0.649 T 0375 x 10 -9 cm 3 s _1 [18]. Equation 
(2.2) is evaluated for the average Lorentz force using the two point function of the magnetic 
field (cf., eq. (1.2)). The final expression depends on the magnetic field strength as well as 
on its spectral index, and can be found in [ 11 ]. 

Dissipation of the magnetic field due to decaying MHD turbulence rests on the fact that 
turbulence is no longer suppressed in the plasma after the decoupling epoch. On scales below 
the magnetic Jeans scale, the magnetic energy on large scales is transferred to small scales and 
dissipates. The dissipation rate of this highly non linear process is estimated by numerical 
simulations of MHD turbulence in flat space using the fact that there exists an appropriate 
rescaling of variables in flat space to match those in an expanding, flat Friedmann-Robertson- 
Walker background [19-24]. In the matter-dominated era, the estimated dissipation rate of 
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a non-helical magnetic field is given by [10] 


decay 


Bq 3 m 


In 


l + Y 


8 vr 2 


In (1 + f: I + In 


1 +z 


m+1 


H(t)( 1 


( 2 . 3 ) 


where Bq is the present-day Held value assuming a flux freezing; m is related to the magnetic 
spectral index as m = • i d j s the physical decay time scale for turbulence given 

by td/U = ( kj/kd ) 5 ~ 14.8 (Bq/1 nG) -1 ^/! Mpc - ) -1 with the magnetic Jeans 

2 2 "fl+ 3 

wavenumber of kj ~ 14.8 n s +5 {Bq/1 nG) n s + 5 {kd/1 Mpc -1 ) n s+ 5 Mpc -1 [10]; and z % < Zd ec 
and ti > tdec are the redshift and time at which dissipation of the magnetic field due to 
decaying MHD turbulence becomes important. 

Ref. [10] calculated the evolution of the matter temperature, T e , and the ionization frac¬ 
tion, x e , with the effect of damping of magnetic fields in the post decoupling era. In ref. [11], 
we revisited these calculations using a modified version of Reef ast++ [25-31]. We found that 
ambipolar diffusion contributes at z < 100, and the effect decreases as hb decreases. The 
amplitude is proportional to Bq. On the other hand, the effect of decaying MHD turbulence 
increases as ng decreases, is proportional to Bq, and dominates over ambipolar diffusion up 
to intermediate redshifts. 

In figure 1 we show the evolution of T e , x e , the visibility function of CMB ( g ), and the 
optical depth to scattering of CMB (r), as a function of redshift. For each choice of parame¬ 
ters, we show the case including the standard implementation of instantaneous reionization 
(denoted as “s r”) and not including it (denoted as “no s r”). We use the CLASS code [32-35] 
to compute these quantities. The standard implementation of CLASS is the same as in the 
CAMB code [36, 37], which uses a tanh fitting formula for the ionization fraction centered at 
the reionization redshift. The case not including the standard reionization is shown only for 
comparison. Dissipation of the decaying turbulence heats the IGM since shortly after the 
decoupling epoch, delaying recombination of hydrogen atoms in the IGM. At low redshifts, 
^ < 20 , the temperature of the completely reionized universe falls as z decreases because 
ambipolar diffusion cannot act on ionized plasma. In a complete treatment of the late time 
thermal and ionization history one should also include contributions from stellar feedback 
which ionizes and can efficiently heat the IGM up to temperatures of 10 4 K at * < 10. Since 
our main conclusions are not based on the temperature but on the ionization fraction, the 
error in the temperature evolution at z < 20 has no consequences. Decay of MHD turbulence 
increases the total optical depth to scattering of CMB shortly after the decoupling epoch, and 
this changes the power spectra of the temperature and polarization anisotropies, as discussed 
below. 

We use the CLASS code with its standard implementation of instantaneous reionization 
at late times to calculate the power spectra, and show the effects of dissipation of fields in 
figure 2 for a magnetic field withBo = 3 nG and spectral indices of tib = —2.9, —2.5, and 
— 1.5. It is well known that the optical depth to scattering of CMB suppresses the temperature 
and polarization power spectra by e~ 2rtot at multipoles higher than that corresponding to 
the horizon size at the epoch of the scattering. In addition to the suppression due to the 
standard late-time reionization, the optical depth due to dissipation of fields, At, creates an 
additional suppression of the power spectra by e -2Ar (see figure 2 ). Also, as the peak of the 
visibility function of CMB shifts slightly toward a lower redshift than the standard decoupling 
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Figure 1: Evolution of the matter temperature (top left), the ionization fraction (top right), 
the visibility function of CMB (bottom left), and the total optical depth to scattering of 
CMB (bottom right), with and without dissipation of magnetic fields, and with and without 
the standard late-time reionization (indicated by “s r” in the legend) of the universe with 
Treio = 0.0925. For better clarity the legend has been included only in two of the figures, 
though it applies to all four figures. To see the effect of dissipation of fields in the standard 
late-time reionization model, compare the blue, red, cyan and pink lines. The cosmological 
parameters are the best fit values of Planck 2013 data plus the low-£ polarization data of 
WMAP 9 [3]. The field strength is Bq = 3 nG with spectral indices of ub = —1.5, —2.5, and 
—2.9. At z > 20, dissipation is dominated by the decaying MHD turbulence. The inset in 
the bottom-right panel shows an increase in the total optical depth due to decaying MHD 
turbulence at z > 200. 


redshift, the angular diameter distance to the decoupling epoch is slightly reduced, shifting 
the locations of the acoustic peaks of the temperature and polarization power spectra to 
lower multipoles. This effect is more pronounced for polarization, and is clearly visible in the 
top right panel of figure 2. 

Not only does an additional optical depth suppress the power spectra at multipoles 
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Figure 2: Temperature and polarization power spectra of CMB, with and without dis¬ 
sipation of magnetic fields, assuming standard late-time reionization of the universe with 
T reio = 0.0925. The magnetic field parameters are Bq = 3 nG for three different values of its 
spectral index ns- The insets show the difference between each choice of spectral index and 
the case without the magnetic held. The additional optical depth to scattering of CMB due 
to dissipation of the fields suppresses the power spectra at high multipoles, while regenerating 
the polarization power spectra in 100 < i < 400. Also, the shift of the visibility function 
to lower redshifts shifts the locations of the acoustic peaks of the power spectra to lower 
multipoles. The cosmological parameters are the best fit values of Planck 2013 data plus the 
Planck 13 likelihood of the low-^ polarization data of WMAP 9 [3]. We also show the Planck 
2013 temperature data [13]. 


higher than that corresponding to the horizon size at the epoch of the scattering, but also 
creates an additional polarization at multipoles corresponding to the horizon size at the epoch 
of the scattering. The amplitudes of the first and second peaks in the F-mode polarization 
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Figure 3: Temperature power spectra of CMB, with and without dissipation of magnetic 
fields, assuming the standard late-time reionization of the universe with r re j 0 = 0.0925. (Left) 
Contributions from the Sachs-Wolfe term (SW), the total integrated Sachs-Wolfe effect (ISW) 
including early and late time contributions, the Doppler term (Dop) and the polarization term 
(Pol). The black lines show the cases with no field, whereas the other lines show Bq = 3 nG 
and riB = —1.5. (Right) The Doppler term for B = 3nG and ns = —1.5, —2.5, and —2.9. 


power spectrum at l « 100 and 400 are enhanced due to regeneration of polarization by 
additional scattering of CMB between the decoupling and the reionization epochs. This 
can be seen in the top right panel of figure 2. While the polarization power spectrum at 
higher multipoles is significantly influenced by the change in the optical depth, that at lower 
multipoles is not. The “reionization bump” in the polarization power spectrum at l < 10 is 
not affected by dissipation of the fields, as ionization at z < 20 is totally dominated by that 
of the late-time reionization. Similar features are observed in the temperature-polarization 
cross power spectrum (cf. lower panel of figure 2). This phenomenology is similar to the 
effect of heating due to annihilation of dark matter particles [3, 38-40]. 

The top left panel of figure 2 shows that the low-f? temperature power spectrum increases 
with increasing magnetic spectral index. We find that this is due to an enhanced Doppler 
term (in Newtonian gauge). For our choice of parameters this is most noticeable for the 
largest spectral index, ns = —1.5. In the left panel of figure 3, we show the contributions 
from the Sachs-Wolfe term, the integrated Sachs-Wolfe term, the Doppler term, and the 
polarization term to the total angular power spectrum. The right panel shows only the 
Doppler term, which increases as ub increases at l < 10. This explains what we find in the 
top left panel of figure 2. 

3 Parameter estimation 

To estimate the magnetic field parameters along with the cosmological parameters, we use 
Markov chains Monte Carlo as implemented in montepython [41]. Specifically, we fix the 
magnetic spectral index, ns, and vary the field strength, Bq, along with the six standard 
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ACDM parameters, i.e., the Hubble constant, Hq = 100 h km/s/Mpc; the physical baryon 
density, u>b = f \h 2 -, the physical cold dark matter density, uj c dm = ^cdmh 2 ', the amplitude of 
the scalar power spectrum at k = 0.05 Mpc^ 1 , A s ; the tilt of the power spectrum, n s ; and 
the optical depth from the late-time reionization, r re j 0 . 

We use the Planck 2013 temperature data [13, 14] (including the lensing data [15]) 
and the likelihood of the WMAP 9-year polarization data [16, 17] derived by the Planck 
collaboration [14]. The lensing data are essential for breaking degeneracy between A s and 
Bq. This degeneracy arises because the optical depth constraint from the temperature data 
is degenerate with A s . (The temperature data effectively measure A s e~ 2Ttot .) In our study 
the key element is that the magnetic held dissipation significantly changes the optical depth 
close to decoupling and so changes the total optical depth to decoupling but does not change 
the optical depth to reionization. The former is determined by the temperature anisotropy 
data and the latter by the polarization data. The combination of the two data sets together 
with the lensing data allows to constrain the magnetic held parameters. 

We summarize the constraints on Ho and the cosmological parameters in table 1, and 
show the marginalized posterior distributions of Bq in hgure 4. We chose hb = —2.9, —2.5, 
and —1.5 as representative cases. We find no evidence for dissipation of the magnetic helds in 
the data. The 95% CL upper bounds on the held strength are Bq < 0.63, 0.39, and 0.18 nG 
for riB = —2.9, —2.5, and —1.5, respectively. 2 

These limits are significantly stronger than those from the magnetic scalar- and vector¬ 
mode contributions to the angular power spectra of the CMB. For example, the Planck 
collaboration used the Planck 2013 temperature data, the WMAP 9-year polarization data, 
and the high-£ temperature data of ACT and SPT to find Him pc < 3.4 nG (95% CL) for 
a smoothing scale of 1 Mpc from the magnetic scalar and vector modes, while ignoring the 
effects on the ionization history of the universe. Our limits on Bq, which is smoothed over 
the magnetic damping scale at decoupling, become stronger if we scale them to a Mpc scale, 
as changing the smoothing scale, X s , to a different smoothing scale, A*, leads to a rescaling 

n B "I” ^ 

of the magnetic field amplitude by a factor of (A s /A*) 2 . 

The only limit in the literature that is stronger than ours is Ho < 0.01 nG for a smoothing 
scale of 0.1 kpc [42], which was derived from the change in the ionization history before the 
decoupling epoch, due to the effect of clumping in the baryon density perturbations induced 
by a primordial magnetic field. Rescaling our limits to 0.1 kpc makes their limit much 
stronger than ours. 

In figure 5 the two-dimensional marginalized posterior probability density distributions 
(Ho versus the cosmological parameters) are shown, for ns = —2.9, ns = —2.5, and ns = 
— 1.5. Most of the cosmological parameters are determined independently of Ho; however, 
the scalar spectral tilt, n s , is weakly correlated with Bq because Ho suppresses the power 
spectrum at high multipoles by e~ 2Ar , which can be partly compensated by increasing n s . 

4 Conclusions 

Dissipation of the magnetic fields in the post-decoupling era heats the IGM and delays 
recombination of hydrogen atoms. This effect can be detected as the extra optical depth 
to scattering of CMB photons. A qualitatively new effect is that it affects only the total 
optical depth to decoupling, which is determined by the temperature data together with the 

2 The corresponding smoothing scales are much smaller than 1 Mpc. The largest one is of order 13 kpc for 
Bq = 0.63 nG, and 4 kpc for Bo = 0.18 nG (cf. equation (1.1)). 
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Figure 4: Marginalized posterior distributions of the magnetic field strength, Bq (in units of 
nG), for the magnetic spectral index of ng = —2.9 (left), n B = —2.5 ( middle ) and ng = —1.5 
(right). 
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Table 1: Best-fit values and 68 % confidence limits on the present-day magnetic field strength, 
Bq (in units of nG), smoothed over k^dec given in equation (1.1), and the standard ACDM 
cosmological parameters. 


lensing data, but not the optical depth to reionization inferred from the low- 1 ! polarization 
data (£ < 10). Using the 2013 Planck data including the CMB lensing data, together with the 
likelihood of the low-1 polarization data from WMAP as derived by the Planck collaboration, 
we find no evidence for the effect of dissipation of (non-helical) magnetic fields in the CMB 
data. The 95% CL upper bounds are Bq < 0.63, 0.39, and 0.18 nG for ng = —2.9, —2.5, 
and —1.5, respectively. These limits are stronger than the previous limits that did not use 
the effect of fields on ionization history of the universe. 

Inverse Compton scattering of CMB photons by hot electrons in the IGM heated by 
dissipating magnetic fields in the post decoupling era leads to a y-type spectral distortion 
[43]. For negative magnetic spectral indices, the Compton y parameter is well approximated 
by (see footnote 1 for the update from ref. [ 11 ]) 

/ R \ 1-7263 

y(n B , B 0 ) = 1.2194 x 10 " 5 \-±j (-ng ) 0 ' 3602 

/ d \ 1-7260 

-1.2155 X 1(T 5 ( ig J (_ nB ) 0 . 3619 e 9 . 3978 xl 0 -%-n B )i 0 . 984 2 ^ 
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Figure 5: 68 % and 95% confidence regions of the field strength, Bq (in units of nG), versus 
the cosmological parameters of the ACDM model. The magnetic spectral indices are ns = 
—2.9 (upper panel), ub = —2.5 (middle panel), and ng = —1.5 (lower panel). 


For the upper limits on the magnetic field amplitudes derived in this paper, the resulting y- 
type distortions are y < 10 -9 , 4x 10 -9 , and 10 -9 for ns = —2.9, —2.5, and —1.5, respectively. 
These limits are well below the contribution from the thermal Sunyaev-Zel’dovich effect in 
galaxy clusters and groups, y « 10~ 6 [44], and thus are negligible. 

We expect that this new method would yield much improved limits, once the full Planck 
data sets including high-7 polarization are used, as they provide significantly better measure¬ 
ment of the CMB lensing (which fixes A s thus breaking degeneracy between A s and e -2AT ), 
and potentially measure the extra polarization at 100 < l < 400 generated by dissipation of 
decaying MHD turbulence. 
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